//  Table A51	  


clear all
clear matrix
set more off

* Data ready for analysis:

use "~/Dropbox/Replication_MVC/Datasets/datasets_analysis/panel_excombatientes.dta", clear 

* Define mean value of economic shock for the group 
sort Exestructura
by Exestructura: egen mean_gold_production_weakties=mean(oro_sin_titulo)
gen mean_gold_shock_weakties=mean_gold_production_weakties*log_gold_price

sum mean_gold_shock_weakties
generate mean_gold_shock_weakties_r = (mean_gold_shock_weakties-r(mean))/r(sd)
 
 
eststo clear
set more off
 
*All captures

reghdfe diffcaptures gold_shock_r mean_gold_shock_weakties_r, absorb(Exestructura year origmun) vce(cluster Exestructura)
estadd local ALLFE = "$\checkmark$"
estadd local TT = "$\checkmark$"
summ diffcaptures if e(sample) ==1
estadd local meanDep = trim("`: display %12.4fc `r(mean)''")
estadd local SDDep = trim("`: display %12.4fc `r(sd)''")
eststo

reghdfe diffcaptures gold_shock_r mean_gold_shock_weakties_r c.pobreza##c.time_t c.pobl_tot##c.time_t c.disbogota##c.time_t c.paved_roads##c.time_t, absorb(Exestructura year origmun) vce(cluster Exestructura)
estadd local ALLFE = "$\checkmark$"
estadd local TT = "$\checkmark$"
estadd local MT = "$\checkmark$"
estadd local BC = "$\checkmark$"
summ diffcaptures if e(sample) ==1
estadd local meanDep = trim("`: display %12.4fc `r(mean)''")
estadd local SDDep = trim("`: display %12.4fc `r(sd)''")
eststo

reghdfe diffcaptures gold_shock_r mean_gold_shock_weakties_r age female black indigenous c.pobreza##c.time_t c.pobl_tot##c.time_t c.disbogota##c.time_t c.paved_roads##c.time_t, absorb(Exestructura year origmun) vce(cluster Exestructura)
estadd local ALLFE = "$\checkmark$"
estadd local TT = "$\checkmark$"
estadd local MT = "$\checkmark$"
estadd local BC = "$\checkmark$"
estadd local CO = "$\checkmark$"
summ diffcaptures if e(sample) ==1
estadd local meanDep = trim("`: display %12.4fc `r(mean)''")
estadd local SDDep = trim("`: display %12.4fc `r(sd)''")
eststo

* Red-handed Captures
reghdfe diffcaptures_f gold_shock_r mean_gold_shock_weakties_r, absorb(Exestructura year origmun) vce(cluster Exestructura)
estadd local ALLFE = "$\checkmark$"
estadd local TT = "$\checkmark$"
summ diffcaptures if e(sample) ==1
estadd local meanDep = trim("`: display %12.4fc `r(mean)''")
estadd local SDDep = trim("`: display %12.4fc `r(sd)''")
eststo

reghdfe diffcaptures_f gold_shock_r mean_gold_shock_weakties_r c.pobreza##c.time_t c.pobl_tot##c.time_t c.disbogota##c.time_t c.paved_roads##c.time_t, absorb(Exestructura year origmun) vce(cluster Exestructura)
estadd local ALLFE = "$\checkmark$"
estadd local TT = "$\checkmark$"
estadd local MT = "$\checkmark$"
estadd local BC = "$\checkmark$"
summ diffcaptures if e(sample) ==1
estadd local meanDep = trim("`: display %12.4fc `r(mean)''")
estadd local SDDep = trim("`: display %12.4fc `r(sd)''")
eststo

reghdfe diffcaptures_f gold_shock_r mean_gold_shock_weakties_r age female black indigenous c.pobreza##c.time_t c.pobl_tot##c.time_t c.disbogota##c.time_t c.paved_roads##c.time_t, absorb(Exestructura year origmun) vce(cluster Exestructura)
estadd local ALLFE = "$\checkmark$"
estadd local TT = "$\checkmark$"
estadd local MT = "$\checkmark$"
estadd local BC = "$\checkmark$"
estadd local CO = "$\checkmark$"
summ diffcaptures if e(sample) ==1
estadd local meanDep = trim("`: display %12.4fc `r(mean)''")
estadd local SDDep = trim("`: display %12.4fc `r(sd)''")
eststo

esttab est* using "~/Dropbox/Replication_MVC/Paper_tables/TableA43_panelA.tex", se ///
	   star(* 0.10 ** 0.05 *** 0.01 ) ///
	   stats(meanDep SDDep N ALLFE TT MT BC CO, fmt(%12.4fc %12.4fc  %12.0fc %12.0fc %12.0fc %12.0fc  %12.0fc %12.0fc) ///
	 label("\small Mean of Outcome" "\small S.D of Outcome" "\small Observations" "\small Municipality, Year, and Group FE" "\small Time Trends" "\small Region $\times$ Year"  ///
	 "\small Municipality Characteristics TT" "\small Individual Covariates" )) ///
		varlabels(gold_shock_r "Economic Shock" mean_gold_shock_weakties_r "Average Shock (Weak)" ) ///
	   ti("Weak Ties Full Results") ///	  
	   replace
	
	
eststo clear
set more off
 
*All captures
reghdfe diffcaptures gold_shock_r mean_gold_shock_weakties_r, absorb(Exestructura year origmun) vce(cluster Exestructura)
nlcom _b[mean_gold_shock_weakties_r]/(_b[gold_shock_r] + _b[mean_gold_shock_weakties_r]), post
eststo 
estimates store All1

reghdfe diffcaptures gold_shock_r mean_gold_shock_weakties_r c.pobreza##c.time_t c.pobl_tot##c.time_t c.disbogota##c.time_t c.paved_roads##c.time_t, absorb(Exestructura year origmun) vce(cluster Exestructura)
nlcom _b[mean_gold_shock_weakties_r]/(_b[gold_shock_r] + _b[mean_gold_shock_weakties_r]), post
eststo 
estimates store All2

reghdfe diffcaptures gold_shock_r mean_gold_shock_weakties_r age female black indigenous c.pobreza##c.time_t c.pobl_tot##c.time_t c.disbogota##c.time_t c.paved_roads##c.time_t, absorb(Exestructura year origmun) vce(cluster Exestructura)
nlcom _b[mean_gold_shock_weakties_r]/(_b[gold_shock_r] + _b[mean_gold_shock_weakties_r]), post
eststo 
estimates store All3

* Red-handed Captures
reghdfe diffcaptures_f gold_shock_r mean_gold_shock_weakties_r, absorb(Exestructura year origmun) vce(cluster Exestructura)
nlcom _b[mean_gold_shock_weakties_r]/(_b[gold_shock_r] + _b[mean_gold_shock_weakties_r]), post
eststo 
estimates store RH1

reghdfe diffcaptures_f gold_shock_r mean_gold_shock_weakties_r c.pobreza##c.time_t c.pobl_tot##c.time_t c.disbogota##c.time_t c.paved_roads##c.time_t, absorb(Exestructura year origmun) vce(cluster Exestructura)
nlcom _b[mean_gold_shock_weakties_r]/(_b[gold_shock_r] + _b[mean_gold_shock_weakties_r]), post
eststo 
estimates store RH2

reghdfe diffcaptures_f gold_shock_r mean_gold_shock_weakties_r age female black indigenous c.pobreza##c.time_t c.pobl_tot##c.time_t c.disbogota##c.time_t c.paved_roads##c.time_t, absorb(Exestructura year origmun) vce(cluster Exestructura)
nlcom _b[mean_gold_shock_weakties_r]/(_b[gold_shock_r] + _b[mean_gold_shock_weakties_r]), post
eststo 
estimates store RH3

esttab est* using "~/Dropbox/Replication_MVC/Paper_tables/TableA43_panelB.tex", se ///
	   star(* 0.10 ** 0.05 *** 0.01 ) ///
	   stats(meanDep N, fmt(%12.4fc %12.0fc)) ///
	   keep(_nl_1 ) ///
	   varlabels( _nl_1 "Peer Effect") ///
	   ti("Economic Shock and Peer Effects") ///	  
	   replace	   

  
		
		 
		
		 
		 


		 
